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Abstract 

The CP-violating phases associated with Majorana neutrinos can give rise to CP-violating 
effects even in processes that conserve total Lepton number, such as /i ^ eee, fie ee and 
others. After explaining the reasons that make this happen, we consider the calculation 
of the rates for the process of the form £ah ia^c and its conjugate iah — > ia^c, where 
£a,(b,^c denote charged leptons of different flavors. In the context of the Standard Model 
with Majorana neutrinos, we show that the difference in the rates depends on such phases. 
Our calculations illustrate in detail the mechanics that operate behind the scene, and set 
the stage for carrying out the analogous, more complicated (as we explain), calculations for 
other processes such as fi s- eee and its conjugate. 



1 Introduction 

Sometime ago we introduced a prescription for identifying a minimal set of CP- violating 
parameters in the lepton sector that are invariant under a rephasing of the fermion fields. The 
prescription holds for any number of fermion generations and, more importantly, accommodate 
the case in which the neutrinos are Majorana particles. Recently 0, based on that work we 
analyzed the dependence of the squared amplitudes on the rephasing-invariants of the lepton 
sector with Majorana neutrinos, for various lepton- violating or lepton-conserving processes, 
giving special attention to the dependence on the extra CP-violating parameters that are due 
to the Majorana nature of the neutrinos ^• 

It was widely believed that these extra parameters appear only in lepton number violating 
processes [|], and there is a lot of discussion in the literature about their possible observable 
effects@, IT I, In], However, in our recent work |]||, we showed that they can appear 



in lepton-number conserving processes as well. The true condition for the occurrence of these 
parameters in a given process seems to be the violation of lepton number on any fermion line 
in the corresponding diagrams, and not necessarily that total lepton number be violated by the 
process as a whole. In processes that conserve the total lepton number, there are in general 
diagrams in which the individual fermion lines change the lepton number, but do so in such 
a way that the changes between different lines cancel in the overall diagram. The interference 



terms produced by such diagrams contain the extra CP- violating parameters that exist due to 
the Majorana nature of the neutrinos. 

However, this analysis was carried out by considering various generic physical processes, 
classified according to whether they conserve total lepton number, or by how many units they 
violate it, and then finding their generic dependence on the rephasing- invariant CP- violating 
parameters. The question of what is the mechanics that operates in a specific process to give 
rise to such effects was not considered there. This question is important when we attempt to 
consider the difference in the rates for, for example, — > eee vs /x — > eee, due to the CP violating 
phases. The issue here is that, unless the final-state interactions are taken into account, the 
calculation of the two rates will be equal (by the CPT theorem) in spite of the fact that CP 
may be violated. Thus, while the arguments and analysis of Ref. |^ are indicative, a specific 
calculation of the effects depends in general on the kinematical and dynamical aspects of the 
particular process considered. 

In order to fill this gap, we consider in the present paper the calculation of the difference of 
the rates for some leptonic processes and their conjugate ones, in the context of the Standard 
Model with Majorana neutrinos. The processes that we consider have the virtue that they have a 
two-body final state, which makes it simpler to take into account the final-state interactions, yet 
they contain all the ingredients to understand and illustrate the issues that we have mentioned. 
In addition, some of the formulas that we will present on the way, are also required ingredients 
in the corresponding calculations for other processes such as — > eee. 

In Section |2| we make some general remarks about the type of process we consider, explain 
why the effect can be seen in some of the processes and not in others, mention the need to 
consider the effect of the final state interactions, and set the stage for the calculations that 
make up the rest of the paper. The effect of the final state interactions, which show up as an 
absorptive term in the total amplitude, is calculated in Section ^, for both the direct process and 
its conjugate. Based on those results, in Section § we compute the difference of the differential 
rates for the process and its conjugate, and thus we are able to show explicitly that it is given in 
terms of the rephasing- invariant CP- violating parameters for Majorana neutrinos. Our outlook 
and conclusions are given in Section ^. Four appendices contain some of the details of various 
stages of the calculation, including a Fierz transformation formula used, the Cutkosky rules 
employed to determine the absorptive term due to the final state interactions, and the phase 
space integrals over the intermediate states required to implement the Cutkosky formula. 



2 General Remarks 

We want to consider processes of the generic form 

44 ^ 44 (2.1) 

where a, h, c, d take values from anyone of the lepton flavors e, ;U, r, with the condition that 
they are not all equal. Such processes conserve total lepton number but in general violate the 
individual lepton flavors. By the usual substitution (crossing) rules, our considerations also 
apply to those processes that are related to these by crossing. 

We classify the processes into two groups, depending on whether or not the diagrams that 
contain the one-loop photon or Z vertex functions also contribute. We denote by Group I the 
set of processes for which the only diagrams that contribute are the box diagrams, as shown in 
Fig. |l[ The requirement for the photon and Z vertex diagrams to be absent is that 

4/4 and 4y^4 (2.2) 



2 




Figure 1: The physical diagrams for the general process iaib — *■ ^c^d- The unlabeled vector boson 
lines represent the W vector boson. In addition to these diagrams, there are other diagrams in 
which any one of the W bosons are replaced by their corresponding unphysical Higgs particle, 
plus the exchange diagrams in which ic and interchanged. Note that diagram (B) contributes 
only if the neutrinos are Majorana particles. 



and 

£d^£a and 4/4 (2.3) 

Because each of the four indices a,b,c,d can only take three possible values (e, /i,r), it is clear 
that the conditions in (|2.2|) and (p.3[) can be satisfied simultaneously only if either 



4 = 4, (2.4) 

or if 



4 = 4. (2.5) 

Without any loss in generality, we can fix one or the other of the two possibilities, and thus 
define the processes in Group I as those of the general form 

44^44 (4/4/4). (2.6) 

On the other hand, the processes in Group II are those for which one (or more) of the 
conditions in and ( |2.3| ) is not satisfied. Again, without loss of generality, they can be 

represented by the general form of Eq. (|2.lJ ) with 4 = 4, i.e., these processes have the generic 
form 



44 ^ 44 , (2.7) 

with c / 6 being the only restriction. In this case, in addition to the diagrams shown in Fig. |^, 
the diagrams shown schematically in Fig. ^ must be included. 



2.1 No-go result for Group I 

As already stated, the diagrams that contribute to the amplitude for these processes are the 
box diagrams shown in Fig. |l|. To the leading order (1/M^), only the diagrams involving the 
W exchange are important. Let us consider 44 ~^ ^-c^c- The physical amplitude is of the form 

A[£a{m{p) ^ 4(pi)4(P2)) = Ma(4(A:)4(p) ^ 4(pi)4(P2)) 

+ Mb(4(A;)4(p) ^ 4(pi)4(P2)) (2.8) 
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Figure 2: Representation of the collection of diagrams for the processes in Group II, that involve 
the one-loop 7^fe^c or Zli,lc vertex functions. 

where Mx{ia{k)ib{p) ~^ 4(pi)4(P2)) is the contribution from each diagram. Their calculation 
is straightforward, and some of the details are provided in Appendix |^ The results for the 
present case are summarized by the formulas 

MB{iaik)£bip) ^ 4(P1)4(P2)) = X^B^"\ucip2h''Lub{p)][uc{pih^,Lua{k)] , (2.9) 



where 



.4 



^5^"'^'^ " E(^a*^^-^-K.,)(4/j:^-^) + (c ^ d) , 

W ij 

The functions /J^'^ and /^^-''^ appearing in these equations are given by 

JA - 




/F = -^<!-f^--.+(^-J)^ (2.11) 



where 



(2.12) 



Further, the Vai are the elements of the lepton mixing matrix, and the Ki are the phases defined 
by the Majorana condition 

r/f = i^^.. (2.13) 

Although there are several equivalent ways to write the results for the amplitudes M^^^, they 
can be brought to this form by suitable Fierz transformations. 

The upshot of this is that when Eq. ( |2.9| ) is substituted in Eq. ( |2.8| ) , the effective couplings 
Xa,b appear in combination as a common overall factor of the total amplitude. Since there is no 
interference term, the rate for the process and its conjugate is the same, and there is no observable 
CP violating effect in this type of process. We have considered explicitly the amplitude for 
44 (^c^c, but similar arguments hold for the inverse 4^c — > 44; whose amplitude is simply 
the complex conjugate, and other related processes such as la ib^^c^c- 
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While we have summarized the result of the actual calculation, a little thought reveals 
what is going on. Because we are calculating to leading order in 1/M^, the dominant terms 
come from the T^-exchange diagrams, as we have already mentioned. The chiral nature of the 
W interactions dictate that, to leading order, only the left-handed components of the external 
fermion fields enter in the amplitude. The most economical way to express this fact is by writing 
down the effective Lagrangian for this process which, by the above argument, can only be of the 
form 

■^^^^ = ^[icl^Lh][ecl^.Lia] + h.c. (2.14) 

In fact, the results given in Eqs. ( ^.8| ) and (|2.9| ) can be represented by this Lagrangian, with the 
identification A = A^^'^'^'' + \^^'^'^\ Thus, to this order the effective Lagrangian actually consists 
of only one term, and therefore the rates for the process and its conjugate are equal. 



2.2 Evasion for Group II 

By the same argument, it is now easy to see how the processes in Group II differ. We consider 
specifically those with 7^ ^j. The diagrams for a process of the generic form given in Eq. ( |2.7| ), 
include the diagrams that involve the ^ibf-c and Zi\,ic one-loop vertex functions [^]. Instead of 
Eq. (|2.8D, the physical amplitude in this case is of the form 



A(ia(k)h{v) ^ U^Yc{v)) = (Af '^^^ +Ai"'"^))[n,(A:')7^^^a(A;)][^c(p')7M^^6(p)] 



where 



with 



(2.15) 



Ar^-^EK.K./?\ (2.16) 



f^P=rk\o^rk. (2.17) 
while X and Y are the neutral-current couplings of the lepton 

X = +sm^ew 

y = i. (2.18) 

We mention the following. In the formula quoted in Eq. ( |2.16| ) we have neglected the other 
terms of order 1/Myy that do not contain the logarithmic factor logrj. In addition, none of 
the terms that arise from the diagram that involve the photon vertex function contain that 
logarithmic factor, and therefore we have omitted altogether that contribution in Eq. ( 2.15| ). 



Thus, by the same argument that led us to write Eq. (2.14), this amplitude corresponds to an 
effective Lagrangian 

^(H.+z) ^ (A(,"''"^) + A(,"''^^))[47^i^4][47^L4] 

+xP[ia7''{X + Yj^KWcl^m + h.c. . (2.19) 
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This in turn can be written in the more compact form 



^(W+z) ^ + y'75)4][47/.^4] + h.c. , (2.20) 

with 

Y' = -{Xf"'"'^ + X^^'"'^) + X^Y . (2.21) 

In particular, X' and Y' are complex quantities, and depend on the Majorana phases through 
the Xb term. 

A kinematic observable that depends on the interference term lm{X'*Y') will be sensitive 
to the CP-violating Majorana phases. However, the total rate, determined from this effective 
Lagrangian, or equivalently calculated with the amplitude given in Eq. ( 2.15| ), will not depend 



on that interference term. The reason is that, as seen from Eq. ( |2.2[l| ), such term arises from the 
interference between the vector and axial vector parts of the current, and that vanishes after 
summing and averaging over the polarizations and integrating over the phase space. Moreover, 
the rates for any process and its conjugate, determined from Eq. ( 2.20D , are equal. This is 



ultimately due to the fact that, while CP does not hold in Eq. (2.20), CPT does hold and that 
is sufficient to guarantee the equality of the rates. 

When we consider the final state interactions between the outgoing leptons, this is no longer 
true as is well known. While the amplitude for the direct process depends on X' and Y' , the 
amplitude for the conjugate process depends on the complex conjugates of these two quantities. 
However, the final state interactions induce an extra phase that is the same for both the direct 
process and the conjugate. This mismatch between the two sets of phases in both cases leads 
to the inequality of the total rates. In the language of the effective Lagrangian, the effect of the 
final state interactions is to augment Eq. ( 2.20| ) in a way that renders it non-hermitian. Our 



task in the next section is to calculate the effect that we have just outlined. 

We would like to remark that the above statements about the equality or inequality of the 
total rates, as the case may be, apply to the total differential rates as well. The latter quantities 
are defined from the squared amplitude as usual, by summing over the final spins and averaging 
over the initial ones, with the total integrated rates being obtained from them by carrying 
out the only non-trivial integration, over the azimuthal angle. The reason why we can also 
consider the total differential rates is the following. Under a CP or CPT transformation, the 
amplitude is related to the amplitude for the conjugate process, with perhaps the momentum 
and/or spin variables reversed. For the total integrated rates, the reversal of the momentum and 
spin variables is of no consequence since they are being summed over. In the case of the total 
differential rates, the reversal of the spin variables is not relevant either for the same reason. On 
the other hand, for the two body processes that we are considering, the momentum vectors can 
appear only through the scalar variables formed out of the scalar products among them, and 
those variables are unchanged by the simultaneous transformations of the momentum vectors[|. 

^The situation is different if we consider for example, processes with three particles in the final state. In that 
case, the spin-averaged squared amplitude can contain momentum contractions involving the four-dimensional 
antisymmetric tensor, and those change when the momentum vectors are reversed. Our considerations could be 
applied to such cases also, but only if some (non-trivial) angular integrations are made so that those terms do 
not appear in some specially designed differential rates. 
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Figure 3: Schematic representation of the one- loop amphtude for ~^ ^c^d- The diagram 
stands for the cohection of diagrams referred to in Figs. |l] and 




Figure 4: Schematic representation of the correction to include the electromagnetic final state 
interactions in iaib (-c^-d- 

3 Final state interactions 

3.1 General considerations 

We denote by M^^^ta^b ^ 44) the amplitude for £Jb ^ 44 determined from Eq. ( ^ ), 
which we represent schematically in Fig. ^. Analogously, we denote by M^^\ia^b ~^ ^a^c) the 
amplitude for the conjugate process. As we have already mentioned, since the individual terms 
that contribute to M^^^ do not contain an absorptive part, then the rates T{iaib — > ^a4) and 
r(^a^f) ^a^c) calculated with the above amplitudes are equal, as a consequence of CPT. This 
result does not hold if we include the final state interactions between the two leptons, of which 
the dominant one is the electromagnetic interaction. Diagrammatically, the additional terms 
are represented in Fig. ^ The total amplitude is 

M = m(°)+m(^4 (3.1) 

Using the Cutkosky rules, the absorptive part of M^-^^ can be expressed in terms of essentially 
M(o) itself, times some factors. 

In order to write the following formulas in a compact form, let us introduce the following 
notation. We denote the initial and the final states in the process iah — > ^a^c by 

\i) = Sa)h{P, Sb)) 

I/) ^ \ea{k',s'Mp',s',)) . (3.2) 
The 5-matrix element for this process is then written as 

{f\S\i) = (27r)^5(^)(fc + p-A;'-p')M°^(44^44)+iM«(44^44)l , (3.3) 



where M^^\£aib £a^c) is determined from Eq. ( 2.20| ). Our next task is to find the appropriate 
expression for the two-loop amplitude M^^^lah — > £a^c)- 
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In general, it can be decomposed as 



m(i) (44 ^ 44) = Mjf ) (44 ^ 44) + (44 ^ 44) , (3.4) 



where M^^^Jp and m''^^ stand for its dispersive and absorptive parts, respectively. Since, as 

already mentioned, M^^^ does not contain an absorptive part, the contribution from M^^^^ only 

produces an order a correction to M^^^ which we neglect. On the other hand, although m'^^ 
is also a factor of order a smaller than M^^^ , it is the piece we are after since it contributes an 
absorptive part to the full amplitude. To calculate it, we employ the Cutkosky rules [0]. As 
shown in Appendix in our notation for the present case they yield 

^ibs(44 - 44) = ^5:(27r)^<5(^)(gi+g2-A.'-p')^^^^^(^i"^4"^ - 44)m(0)(44 - ^i"^4"^) • 

n 

(3.5) 



Here the symbols ih^ stand for an intermediate 4^c state 

I"-) = |4(gi,cri)4(Q'2,cr2)) , (3.6) 

and the sum over the intermediate states stands for 



The quantity M^''') in Eq. ( |3.5| ) is the electromagnetic scattering amplitude for 4('i'i)Cri) + 
4(^2, 02) £a{k', O + Up', s'c), i.e.. 



zM(^)(^i")4") ^ 44) = (^^^^ [Ua(fc',sl)7'^^Z„(gi,CTi)] [^,(p',s;)7^t.e(92,CT2)] . (3.8) 

Thus finally, the full amplitude is given by 

M(44 ^ 44) = m(o)(44 ^ 44) + M^^{iah ^ 44) , (3.9) 

where M^bs(44 ^ 44) is computed from Eq. (^^, with M(°)(44 ^ 44) determined from 
@)- 

In a similar fashion, for the conjugate process 44 ^a^c, we define 

\i) = \Lik, Sa)ib{P, Sb)) 
I/) = \ia{k',s'Mp',s',)) 

\n) = 14(91,^11)4 (92, fT2)) . (3.10) 

Neglecting again the dispersive part of the amplitude, the corresponding 5'-matrix element is 
given by 

i) = {2Tr)^6^^Hk+p- k' - p') [iM(o)(44 ^ 44) + ^mS^(44 ^ 44)] . (3.11) 

where 

^ibs(44 ^ 44) ='-Y.{2TT)H^''\qi+q2-k'-p')M^^\l^;^nt) ^ 44)mW(44 ^ ' 

n 

(3.12) 
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with 

iMW(^W4") ^44) = ^q'^lyy [va{qi,a^)l^Va{k',s',)] %{q2,a2hMp' ^ ■ (3-13) 

It is useful to note the foUowing. Using the relation between the spinors u{p, s) and v{p, s), 
e-g-, 

v{p,s) =i-f2U*{p,s) , (3.14) 
in a specific convention, together with relations such as 

v'^jf^v = u'jf.u' , (3.15) 

it follows that 

^ Uc) = ^ 44) , (3.16) 

that is, the electromagnetic amplitude for the two processes is the same. 

Our task at hand is to apply these formulas to compute the absorptive part of the amplitudes 
for the direct process and its conjugate, from Eqs. (|3.5|) and ( ^.121) , respectively. 

3.2 The absorptive part m£^(44 ^ 44) 

First of all, from Eq. ( |2.20| ) we write 

mW(44 ^ ^i"^4"0 = [ua{qih^{X' + rj,)uaik)][u,iq2h^,Lub{p)] . (3.17) 
Using this and Eq. (|3.8D , we then obtain 

^ib!(^a4 ^ 44) = ^ / di[tia(A:')7\^i + mahf^iX' + y'75)^^a(A:)] 

X [ucip'hxih + mchf,Lub{p)] , (3.18) 

where the symbol J dL stands for 

^ (27r)324f (27r)324^) ^ ^ ^ (fc' - q,)2 

In order to proceed, we introduce the following definitions for the integrals over the intermediate 
momenta 





j,lL 






jiL 


qifi 




jiL 


q2ti 


/(12) _ 


jiL 





(3.20) 



in terms of which 



^ibs(^a4 - 44) = + + ^4'^ + ^4'^ ) , (3-21) 
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where 



Mi'^ = marric 



M^'^ = mjj^'^\ua{k')j^Yr{X' + Y'j5)ua{k)\ \uc{p'hxlf.Lu,{p) 



Uaik'h^^nX' + Y'^5)Ua{k) Ucip'hxl" IpLu.ip) 



(3.22) 



While the integrals are doable in the general case, the procedure is tedious and the final formulas 
are cumbersome. Therefore, for the moment we proceed as far as possible without using the 
explicit results of their evaluation. 

3.3 The absorptive part M^ll{Ub Uc) 

We carry out the same procedure with the amplitude for the conjugate process. From Eq. (^3o|) , 



(3.23) 



which, using relations such as those give in Eqs. ( 3.15| ) and ( |3.14| ), can be written in the form 



M(0)(44 ^ = \uaiqi)Y{X'* - Y'*j,)uaik)] \uc{q2h^.Rub{p) 



(3.24) 



From Eq. ( ^.12 ), and using Eq. ( |3.16| ), we then obtain 



Mi^sM ^ Uc) = ^ I dL ua{k')^\ii + ma)r{X'* - Y'*^^)Ua{k) 



Ucip'hxih + mc)-f^Rub{p) 



(3.25) 



By comparison, it is immediately seen that the amplitude for this process is obtained from the 
formulas for the direct process by making the substitutions 



(3.26) 



X' - 


X'* 


Y' - 


-> y'* 


75 - 


-75 



4 The difference in the rates 

If we write the total amplitude in the form 

M = M^^\lJb ^ + M^^^iiJb 
M = M^''\Ub ^ Uc) + M^^iUb 

the quantity in which we are interested is the difference 



Ml 



2 Re (MW*(44^44)Mi^)(£ 

2 Re (mW*(44 ^ Uc)Mi^^{Ub ^ U) 



(4.1) 



(4.2) 
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where the angle bracket notation indicate the operation to sum and average over the final and 
initial spins, respectively. In this expression we have made use of the fact that, after that 
operation is made, the terms without the absorptive part cancel out. Using Eq. ( 3.21| ), and 
denoting M^^^ {ia-^b ~^ ^a^c) by simply M^^^ , we can write 



Ml 



R — R . 



where 



(4.3) 



(4.4) 



i=l 



and R is obtained from R by making the substitutions indicated in Eq. ( 3.26| ). We now compute 
the various terms in R. 



4.1 1st term 

Averaging over initial spins and summing over final spins, we obtain 

1 



4 

X Tr 

_ 1 
~ 4 
X Tr 



mamc 



j{0)Tr 



(l^ + mb)-faL{'^' + mchxjf^L 



+7'^^V7^(|X'|2 + |y'|2 + X'*Y'j5 - Y'*X'-f5 
The traces are easily evaluated with the help of the formulas 



Tr7„ 7/37^,75 
Tr 7a 7/377 7575 



where 



Co/375 — gapQyS — 9a^9(}S + OaSgpy 



(4.5) 



(4.6) 



(4.7) 



(4.8) 



When Eq. ( |4.7D is used in (4^), four terms are produced, which we schematically denote as CC, 
Ce, eC and ee, indicating which factor C or e they contain from each of the two traces that 
appear. It is easy to see that the term CC is real, while the terms Ce and eC are zero after 
contracting the corresponding Lorentz indices. Only the term ee has a non-zero imaginary part, 
and a little bit of algebra shows that 



Im (m(°)*m|^^) = 2Amlmll^^h' ■plra{X'*Y') . 



(4.9) 



4.2 2nd term 

For the remaining terms, it is useful to use the identity 



(4.10) 
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where C'^p^^ is defined in Eq. (14.81). It then follows that 



(4.11) 



and 



4 
X Tr 



X Tr 



(4.12) 



It is not difficult to see that all contributions involving the C term from the first parenthesis are 
real. Among the other terms that are not necessarily real, some are zero identically after con- 
tracting the Lorentz indices and there are others that, while not zero identically, are proportional 
to either one of the following factors 



(4.13) 



(2) 

Since the integral Ip is a vector that depends on p' and k', it is proportional to either p'p or k'p. 
Whence all such terms eventually appear contracted in the form 



Apo-T 



kpPaP'rk'x , 



(4.14) 



and in the end yield zero by momentum conservation. In summary, in the first trace that 
appears in Eq. (4.12), only the terms that contain the combination X'*Y' — X'Y'* contribute 
to the difference in the rates, and a little bit of algebra yields 



Im (m(°)*M^^^\ = -8mllm{X'*Y')l''^^p\k' p ■ p' + ppk' ■ p' + p' k' ■ p 



(4.15) 



4.3 3rd term 

Using Eq. (|1^) once more, by straightforward algebraic manipulations we obtain 

1 



-mJ^P X Tr 
4 P 



(|i + mi,)-/pL{^' + nicj-fx^pL 

I + m,)7^(X'* + y'*75)(r + "ia)7a(^' + ^'75 
(^ + maW{X'* + y'*75)(r + ma)lc.{Y' + ^'75) 



(4.16) 



The remaining manipulations and arguments are similar to those that lead to Eq. ( 4.15| ), and 
in this case they yield 



Im l^M^^^M^^^ = 2^'m\ml\^(X'*Y')lf'^pP . 



(4.17) 



4.4 4th term 



It is useful to notice that the expression for Uf^ can be simplified by using the identity in Eq. 
(P^ ), and using then the formulas 



e 



(4.18) 
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Thus we obtain 



M, 



(1) 



and then 



+ li^^\x' + Y')[u,{k'h''Ruaik) 



Tr [{^ + mh)-f^L{^' + m,)7"L] 

+ y') Tr [(^ + mj7^(X'* + y'*75)(r + rriah^R 



(4.19) 



(4.20) 



When we carry out the traces and perform the contractions, the terms that do not have the 
factor of turn out to be real, proportional either to |X' — y'p or to |X' + y'p. The only 
terms that contribute to the rate difference are those that have the factor of m^, and by the 
same manipulations that lead to Eq. ( 4.12[ ) we find 



Im (MW*Mi')) = 8m2lm(X'*y')^^'^ 



pap'f^ + p'"/ + g^Pp . p' 



(4.21) 



4.5 The sum 

Summarizing the formulas that we have obtained, we can now write 

R = -e^ Im {X'*Y') [Zi + Z2 + Z3 + Z4] , 

where 

Zi = 24mlmll^°^p ■ k' , 

Pp{p' ■ k') + k'p{p ■ p) + pp{p ■ k') 



(4.22) 



Z2 = -Smll^^'^P 
Z3 = 24m>2j(i)pP 



(4.23) 



It is now obvious that R, which is to be computed similarly but with the substitution indicated 
in Eq. ( |3.26| ), is given by = —R, and therefore 



|M|' 



2R. 



(4.24) 



Clearly the Zj, and consequently the CP violating effects given by R, vanish in the limit that 
all the charged lepton masses are taken to be zero. We then consider the quantities Zi evaluated 
to the lowest order in the charged lepton masses; i.e., we keep only those terms that contain 
two powers of the charged lepton mass. At this order, Zi and do not contribute. Since Z2 
and Z4 already have an explicit factor m\, we evaluate the other kinematic factors, for massless 
particles. As shown in Appendix in this limit the relevant integrals are given by 



/{2) 
/(12) 



1 r 



BiQp 

{Bo - B2)sg^, + {2Bo - Bq + 52)P^P, 



+ 2Bi{Q^P, - P^Q,) - (3^2 - Bo)Q^Q, 



(4.25) 



13 



where 



P = k' + p , 
Q = k'-p', 
s = P\ 



and 



Thus 



-Bri 



1 



1 tr. 



IGvrs J -I 1 — 



(4.26) 



(4.27) 



Z2 = -Smls [Bop ■ P - Bip ■ Q] 

^4 = imls [2Bop ■P-iBo + B2)p ■ Q] , 

and we finally obtain 

Z = Z2 + Z4 = Amls[2Bi - Bq- B2]p-Q + 0{mjs) 



(4.28) 



(4.29) 



It is reassuring to observe that, while the integrals Bn defined in Eq. ( |4.27| ) are (infrared) 
divergent individually, the combination that appears in Eq. (4.29) is divergent-free, and its 
value is given by 



2Bi - Bq - B2 



8lTS 



Thus, from Eqs. (|]2|) and (|]29|), 

R = 

and finally, from Eq. ( 4.24| ), 

|M|' 



=2 
2^ 



m2lm(X'*y')(p-Q), 



Ml 



-m\ Im {X'*Y'){p ■ Q) . 



(4.30) 



(4.31) 



(4.32) 



On the other hand, the leading term of the amplitude squared is straightforward to calculate 
and yields 



|m(o)|' 



M 



(0)|2 



4|X' + Y'Y{k ■ p){k' ■ p) + 4|X' - Y'Y{k ■ p){k' ■ p) , (4.33) 



which determines the total rate. Taking the massless limit approximation, and using Eq. ( [4.32| ), 
we then have 



1 
f 



dV 



d(cos 0) d{cos 9) 



(m^/s) Im {X'*Y') 



4tt \x' + y'|2 + (i/3)|x' - y'|2 



cos I 



(4.34) 



where 9 is the angle between k and k', in both the direct and the conjugate processes. 
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5 Discussion and Conclusions 



Using Eq. (|2.21| ) and the formulas given in Eqs. ( p.lOD and ( |2.16D , the CP violating quantity 



Im(X'*y') that appears in Eq. (4.34) is proportional to 



{X + Y)Y^ i/J^V?^ [im (%cfc)|Kd' + Im {ta^cjhjck)\VcJ\ 



(5.1) 



where X, Y are the neutral-current couplings defined in Eq. ( |2.18| ) , fjl'^^ and are the kine- 
matic factors defined in Eqs. (2.11) and ( p. 17 ) while the coefficients tbjcfc and Saij are given 
by 



taibj 



Va^V:,K*K, 



(5.2) 



The t coefficients are the rephasing-invariant parameters of the lepton sector that are analogous 
to those introduced for the quark sector in Refs. |15, 16|. These parameters occur in a purely 
lepton-number conserving theory. In fact, if lepton-number is conserved in the theory, the 
diagrams in Fig. |i]B do not exist, so that we can put f^''^ = 0. In this case, Eq. (5.1) clearly 
shows that only the t-invariants appear in the CP-violating part of the amplitude. 

On the other hand, the s coefficients are precisely the rephasing invariants introduced in 
Ref. ||l| to accommodate the Majorana neutrinos. As shown in there, and further studied in 
Ref. the s coefficients form a suitable set of rephasing-invariant parameters for describing 
the CP violating effects due to the Majorana nature of the neutrinos. In fact, the dependence 
of 6 on the product of three s parameters, as indicated in Eq. (|5.l| ), was anticipated in Ref. H] 
[e.g., Eqs. (3.27) and (3.37) of that paper]. 

Thus, the present calculation confirms the expectation that the extra CP violating phases 
that exist for Majorana neutrinos can appear in CP violating observables in processes that 
conserve total lepton number. Although the type of process that we have specifically considered 
(e.g., e + n ^ e + t) is not a realistic one at present, related processes such asr^e + e + fi 
will show the same effect. The corresponding calculations for the latter kind of process is more 
involved than those presented here due to the three-body final state involved. Nevertheless, the 
present calculations, besides serving as a proof of concept, set the stage for considering such 
three-body decay process, and should prove to be technically useful in that context as well. 
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A Calculation of the box diagrams 

We calculate here the diagrams shown in Fig. ||, for arbitrary incoming and outgoing lepton 
flavors. We denote the momentum vectors by k^, with £ = a,b,c,d according to the labels 
assigned in the diagrams. 
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Diagram (|T[A) 

Not counting the exchange diagrams, we have four diagrams altogether. However, those that 
involve either one or two unphysical Higgs in the internal lines are of order and we 

neglect them. Therefore, in terms of the couplings of the charged current jj^^^ = J2 yai^a^fiLui, 



X [VciV*iUc{kc)j^iLiSu^{q)-fxLuaika)] 



This expression can be written in the form 



-ig 



Xp 



{k,-qY-M^ \{ka-qr-M^^ 



(A.l) 



iM 



where 



(abed) _ f 9 \ 



[-^) Y.(^<ijV,*V,^:tK [ud{kdhu7p7^^Lul>{kk)\ [uc{kch''jVLua{ka)\ ,(A.2) 

^5 J 



d^q 



(27r)4 [{k, - qY - Mlr][{ka - qY - " -ka + q? - ml^] 

Neglecting terms 0(1/M^), 



(A.3) 



IpX = -;9pX 



f. 



Hi) 



(A.4) 



where /^*"'^ is given in Eq. in the text. Finally, using the identity 



Ud{kd)lvlXlp.LUh{h) Uc{kc)-/^-f^Y Lua{ka) = 16 [udikdhxLubikb)] u^kch^ LUaika) 



we arrive at 



where 



M^^"^^ = [ud{kd)-i^Lub{kb)][uc{kc)-ip,Lua{ka)] ■ 
The contribution to the physical amplitude, including the exchange term, is given by 



(A.5) 
(A.6) 
(A.7) 
(A.8) 



Using the fact that Ai^^^^^ = —A4ylf"^'^\ which follows from a Fierz transformation, the contri- 



bution to the physical amplitude is then 

MA{iaib^lcid) = Xf'"'^M^^''''^ 

with as defined in Eq. ( 2. 101 ) in the text. 



(A.9) 
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Diagram (|T[B) 

As with the diagrams ||A, there are four diagrams not counting the exchange diagrams, and 
those that involve either one or two unphysical Higgs in the internal lines are of order 
Therefore, to order 



iM 



(abed) 
B 



Ucikc)i-nuLVcj)iKpSuj{ka - - q){ilpRVdj)vd{kd) 



Vb{h){i-i^RV^^iK'^Sy^{q){-i-ixLV*^Ua{ka 



{h + qf - Ml. 



-ig 



{ka-q?-Ml 

(A.IO) 



where we have used the Majorana condition of Eq. ( p.l3| ). Only the mass terms of the neutrino 
propagators contribute due to the L and R factors on opposite sides of the propagator, and we 
obtain 



iM 



(abed) 



Ucikch^YRVdikd) Vb{kb)-fxlpLUa{ka 



where 



Ub 



(2vr)4 [{k, + q^ - M^][{ka - q? - M^^W " mlMK - k^ - q? - ml 



(A.ll) 



(A.12) 



Neglecting terms 0(1/M^), the explicit evaluation of Ib yields 



Ib 



Ti log Ti 



ri + {i^ j) } . 



(A.13) 



Finally, by a Fierz transformation, 

Uc{kch^YRvd{kd) Vb{kb)^x7pLua{ka) = 2 Ud{kd)^^ Luh{kh) Uc{kc)^xLua{ka) , (A.14) 
the details of which are shown in Appendix ^ Thus Eq. ( A.ll| ) reduces to 



M 



(abed) 
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J2 {K!Kfv:^,*VcjVd,) {2g'^)M 



(abed) 
W ' 



(A.15) 



where A4^^^^^ has been defined in Eq. ( A.7 ). As with the previous diagram, the contribution 
to the physical amplitude, including the exchange term, is given by 



Mb (44 - 44) = Mi^'^^ - 



^ {abed) , , (abed) 
^B -^W ' 



(A.16) 



with x^^^^'^^ as defined in Eq. ( p.lOl ). 
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B Fierz transformations 



Fierz transformations concern products of two fermion bilinears. For arbitrary spinors wi, W2, 
and W4, we denote these products collectively as Cj with i = S,V, T, A, P. In other words, 



with 



= (IjTM'f^Mi^' 7^(75, 75) 



(B.l) 
(B.2) 



and r* being the same things with contravariant Lorentz indices. Then we will denote by e- the 
same quantities, but with W2 <-> w^; i.e.. 



e- = [wiT'w4\[wsriW2] . 
The basic Fierz identity is the relation between the two sets of bilinears, {ei} and {e-}, 

j 

where 



(B.3) 
(B.4) 



F 



1 


1 


1 


1 
2 


-1 


1\ 




4 


-2 





-2 


-4 




12 





-2 





12 




-4 


-2 





-2 


4 




1 


-1 


1 
2 


1 


1/ 



(B.5) 



To prove Eq. ( A.14| ) from here, first note that the identity 



(B.6) 



can be used to write the left hand side of Eq. ( A.14| ) as Aes — gt in the notation of Eq. ( |B.1| ), 
with 

W2 = Vjid{kd) 

m = VRbih) 

W4 = ULa{ka) ■ (B.7) 



The relevant Fierz formula for us now is 

Aes — ex '- 



e's + e'p + ^c't + e'v - 



3e's + Se'p - i^e.^ 



= -2{e's + e'p) +e'T + e'y - e'^ ■ 
Substituting the spinors a-d given above, it turns out that e'g prp = 0, while 

- = e'y = ULc(kch^ULa{ka) Vph{kb)l\VRd{kd) 

= ULc{kc)l^ULa{ka) ULd{kd)l\ULb{kb) , 

using Eq. ( p. 15 ) and similar relations in the last step. This gives Eq. ( A.14| ). 



(B.9) 
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C Derivation of Eq. ( pT^l ) 



In order to make the discussion generally applicable and not tied to any particular channel, we 
consider the process labeled as 



L{ka) + ^ Ukc) + Ukd) 



(C.l) 



Thus, this includes processes such as ia^^a and ia'l^b ^a^c- The treatment for the 

crossed processes such as Ia(-b ^a^a is similar, with the appropriate modifications dictated by 
the usual substitution (crossing) rules. The amplitude M^''), determined from the diagrams that 
are schematically represented in Fig. ^, can be expressed in the form 



iM(°) = i X! 9ABlAB{ka, h, kc, kd) Uc{kc)T AUa{ka) Ud{kd)T BUb{h) 



A,B 



(C.2) 



where each qab denotes a product of coupling constants, each Iab is Feynman integral which 
is real, and the Ta are the generalized Dirac-Pauli matrices. In the most general case, the 
amplitude can be brought to this form by making the appropriate Fierz transformations. With 
this in place, the amplitude for M^^-*, determined from the diagrams represented in Fig. ^, is 
given by 



/ T^^i^ 51 9ABlAB{ka, kb, kc + q,kd - q)iD''/{q) 

A,B 

X [uc{kc){-ie^-t^i)iSFc{kc + q)TAUa{ka)] 
X [ud{kd){-iedlu)iSFd{kd - q)TBUb{kb)] . 



(C.3) 



Our task is to determine the absorptive part of M^^\ 

Since the integrals Iab are real, the absorptive part can arise only from the denominators 
of the lepton propagators in Eq. ( C.3 ), which we write in the form 

D = iq'- a)(gO _ a')(g° - b)iq^ - b') , (C.4) 

where 



E 



a 
b 
b' 



kc+q 

kd-q 



ie 



with 



p 



p'^ +m?f . 



Therefore, we rewrite Eq. ( p.3|) as 



d^q 1 



(C.5) 



(C.6) 



(C.7) 
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a' 
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C 








a 


b 



Figure 5: Original path of integration in the plane. 




Figure 6: Deformed path of integration in the plane to avoid the poles when e ^ 0. 



where 



iM^^^ = i ^ gABlAB{ka, h, kc + q,kd- q)iD^/{q) 
A,B 

X [ucikc){-iec^fj,)i{jtc + ^ + mc)TAUa{ka)] 
X [ud{kd){-iedlv)i{h -i + md)TBUb{kb)] 



(C. 



From a mathematical point of view, the absorptive part of M^^^ arises from the fact that, 
as a function of the integrand in Eq. ( |C.7] ) has poles at the points indicated in Eq. ( p.5| ), 
and illustrated schematically in Fig. ^. As e ^ those poles lie on the real axis but, as long 
as the path of integration can be deformed such that it avoids the poles, the resulting integral 
is real and the absorptive part of M^^^ is zero. This is what happens for all the kinematic 
configurations in which none of the poles coincide with another one, as illustrated in Fig. ^. 
However, if the kinematic configuration is such that, in the limit e ^ 0, one of the poles that lie 
above the real axis coincide with one that lies below the real axis, then the path of integration is 
pinched between the two points, and it cannot be deformed to avoid the poles. In this case the 
amplitude will develop an absorptive part. This side of the coin is also illustrated in Fig. if 
we consider the case that b' = a. As a matter of fact, from the condition that Ep^ is a positive 



quantity, it follows from Eq. (C.5) that this is the only possible pinch condition. 
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Figure 7: Deformed path of integration in the plane. 

In order to isolate the contribution to M^^^ in this situation we proceed as follows. By 
Cauchy's theorem, we can deform the original path of integration as shown in Fig. |^ The virtue 
of this is that we can write 



(27r)4 D 



+ 



Ci 



d^q 1 
(27r)4l) 



(C.c 



C2 



By the same argument that we have explained above, it now follows that the integral over the 
path C2 does not produce an absorptive part because the path cannot be pinched (i.e., there 
is no kinematic configuration for which b can become equal to a' or b'). Therefore, the integral 
over C2 contributes only to the dispersive part of the amplitude and we neglect it. The integral 
over Ci on the other hand, can be evaluated by the method of residues, and therefore 



d\ 



(27r)3 D' 



where 



D' = 2E 



kc+q 



id) 



.(c) 



'kd-q 



id) 



E, 



(c) 

kc+q 



E, 



id) 



+ ie 



(C.IO) 



(C.ll) 



Once again, this has a dispersive and an absorptive contribution. Retaining only the latter, 
which is obtained by the substitution l/(x + ie) — > —iirS^x), we finally obtain 



-vr 



d^q 



1 



1 



id) 



E: 



(c) 

kc+q 



Eltq)^M^'' 



(C.12) 



To write this in its final form, we put in the expression for A4^^^ 



qc 
Qd 



kc + q, 
kd - q, 



(C.13) 



and insert the factor 



: (27r f6^'\c 



kc-q) 



I ^^.i'^^f^^'H^d-ka + q). (C.14) 



(2vr) 
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When these are substituted in Eq. ( C.12| ), the integral over d?q can be ehminated with the help 
of the delta functions, and we finally arrive at 



-27 7 (27r)324c^7 (2^)32^, 
with A4^^^ expressed in the form 



i27Tf6{qc + qd-kc-kd)iM^'\ (C.15) 



Using the relation 



A,B 

X [ucikc){-iecj^,)ii4c + mc)TAUa{ka)] 
X [ud{kd){-ied^u)i{id + md)TBUb{kb)] . 



+ rrii) = ue{qe)ue{qi) , 



(C.16) 



(C.17) 



it is easily seen that the result given in Eq. ( C.15| ) is equivalent to the formula quoted in Eq. 

!)• 



D Integrals over intermediate states 



Here we consider the evaluation of the integrals defined in Eq. ( 3.20 ), the results of which are 
quoted in Eq. ( [4.25| ). For the reasons mentioned in the text, we take all the lepton masses to be 
zero. 



The measure dL, defined in Eq. ( 3.19| ), contains the factor 

1 



{k' - qif 



(D.l) 



where ^ is the cosine of the angle between qi and k' , and s has been defined in Eq. ( 4.26 ). Then, 
for any integrand F, we can write 



dL F 



1 



+1 F 



1-e 



(D.2) 



performing as many integrations as the delta function allows us. Further, complementing Eq. 
( 4.26| ) it is convenient to define 



P' = qi+ q2 

Q' = qi- q2- 

Noting that momentum conservation ensures that P' = P^^, it follows that 



/(I) 



^dL {P'^ + Q'^) = ]^{P^J + J^) , 



and similarly 



(D.3) 



(D.4) 



(D.5) 
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where we define a new set of integrals 

J = J dL 

= / dLQ'^^Q'^^...Q'^^^. (D.6) 
For the integral with two indices, we can similarly write 

Therefore, we only have to evaluate the integrals in Eq. (|D.6| ). 
The scalar integral J can be determined immediately. 



where Bn is defined in Eq. (4.27). For the others, notice that 

=0 (D.9) 
for any i = 1, 2, • • • n. For the one-index integral, there is only one such relation 

which implies that 

J^ = bQ^, (D.ll) 

for some invariant b. The invariant can be determined by contracting both sides with Q^, which 
gives 



By explicit computation. 



and using = —s, we arrive at 



6=^. (D.12) 

g^J^ = -2sBi , (D.13) 
J^ = 2BiQ^. (D.14) 



Substituting this in Eq. ( |D.5| ) then yields the formula for ij^^ quoted in Eq. (4.25). An analogous 
formula for I^^^ follows from Eq. (p.4D. 



For the two-index integral J^,y the relation in Eq. ( p.9| ) now dictates the general form 

Jf,^ = KP^^Pu - sg^u) + pQ^Qv , (D.15) 

with the coefficients being easily determined by contracting with g^^'^ and by Q^Q'^. The equa- 
tions obtained this way are 

^s\ + sp= -g^'^J^^. (D.16) 
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The combinations on the right side of these equations are computed exphcitly as 

Q'^Q'^Jf.u = JdL{Q- Q'f = 2s^B2 . (D.17) 
Substituting these in Eq. ( D.16| ) and solving for A and p, 



J„u = (Bo - B2)(P^,Pu - sg^,^) + (3^2 - Bo)Q^Q^ , (D.18) 
which together with Eq. ( D.7 ) yields the formula for iju^^ quoted in Eq. ( [4.25 ). 
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